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XIII. On a new Method of Approximation applicable to Elliptic and Ultra-elliptic Functions. 
By Charles W. Merrifjeld. Communicated by the Rev. H. Moseley, F.R.S. 

Eeceived March 26,— Bead May 24, 1860. 

The difficulty of finding approximate values of elliptic functions of the third kind has 
led me to consider a general method of approximation, which I believe to be new, at 
least in its application to the evaluation of integrals of irrational functions. 

It depends on the known principle that the geometric mean between two quantities 
is also a geometric mean between their arithmetic and harmonic means. If we take 
any two positive quantities, we may approximate to their geometric means as follows : — 
Take the arithmetic and harmonic means of the two quantities, then again take the 
arithmetic and harmonic means of those means, and so on: the successive means will 
approximate with great rapidity to the geometric mean. 

To judge of the convergence of the method, I give, in the first two columns of the 
following little Table, the arithmetic and harmonic means thus derived from the numbers 
1 and 2 (which is the most unfavourable case that need present itself). The third 
column contains the difference of the first two, within which lies the error of either. 

3 4 1 

2 3 6 

17 24 1 



12 17 204 

577 816 1 



408 577 235416 

665857 941664 1 



470832 665857 313506783024 

Either of the fourth pair would thus give the square root of 2, correct to eleven places 
of decimals. 

This method finds its application in evaluating the integral §fx.\/(px.dx, where fa 
and <px are rational functions of x which both increase or decrease regularly, and have 
no singular values, within the limits of integration. If we find successive means, as 
above described, say A^ and H A ,, between fa and fx.(px, then it is clear that, since 

fa\/<px always lies between A* and H,,, so §fa.\/(px.dx must lie between §A w .dx and 

JH^. dx. Now A, and H^ are each the product of fa and a rational function of <px, and 

are therefore themselves rational functions of x. They are therefore always integrable. 

It is not to be denied that the application of this method to the functions which call 
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for its use is cumbrous ; but there is no other way of treating them which is not open 
to the same complaint, for, in reality, they are exceedingly complex of themselves. 

The following are the first three terms of the arithmetic means, with the correspond- 
ing harmonic means written under them, derived from x and y : — 

x-\~y x* + 6xy 4- y* x 4 + 28 ^+7Q^V + 28^ 3 ±J 4 

2 ' 4(a? + y) "' ~8{^ + ~7x*y + 7xy* + f) ' 

2xy 4xy (x -f y) Sxy (x s + 7x % y 4- 7xy® -f y s ) 

oc + y 'x J T<^xy~+'y % " x A +l^'V^^V^2^^Tp* 

It is convenient to use the arithmetic in preference to the harmonic series for inte- 
gration, as the divisor contains one binomial factor less. The third arithmetic mean 
may be resolved into the following form ; — 



xJr V l X V * ( 2 + ^2) an/ 1 (2— \/2)ory 

' WWnt » » « """"* '" ' lll - J ' —1— "" ' ' ' - - ■■ ' -*■ wr .,i. m~ > *m~* m~ ~ ~*-~~ .,■■■■ ~ ~~ ~~- *mm-*, , +., — — • 



Note that 



8 ' 2 w + y l 2 x+{S + 2V2)y ! 2 x+(3-2 \/2)y 



(2+^2) 2 =2(3+2 N /2): 



and that 



3-2 V2 (2- v / 2) 1 



log(3+2 N /2)=:0*76555 13706 75726. 

If a further approximation be thought necessary, it is possible to resolve into partial 
fractions the fourth arithmetic mean ; but if we go beyond this we shall have to solve a 
reciprocal equation of the eighth degree with all its roots real. The fourth arithmetic 
mean will have for its divisor 

the roots of the biquadratic factor, with their signs changed, are 

7+ 4s/2 + i/(80+56>/2) and 7+4^/2- ^(80+56^2): 
their approximate values are 

25-27414 23690 882 

2-23982 88088 434 

0-44646 26921 718 

0-03956 61298 966 

The difference between any (say the nth) pair of means of the series has always {x—yj 
to the power of 2 n for its numerator, and the product of the denominators for its deno- 
minator. The logarithm of the error is therefore always much less than 

2- log (^— y)— (2»— 1) log (x+y)-nlog 2. 
I shall now indicate the mode of applying this method to the general form of an 
elliptic function, J j+^4 ^I=TO^ =^- 

It is obvious that the nearer the ratio x : y is to unity, the less number of terms shall 
we require to obtain a given degree of accuracy. In elliptic functions which involve a 
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radical of the form \/l — <? 2 sin 2 <p, this ratio may always be brought within the limits of 

12 1 

2 and j ; for, if c sin <p should exceed -^=, we may put the radical under one of the forms 

cos (p \/l+(l—c 2 ) tma 2 <p or \/l—<? 2 /\/l4*Yzr2tan 2 (p, and one or other of these new 

radicals will always be less than \/2. I should have remarked that c never exceeds 
unity. 

We may take u as a geometrical mean between v and ^(1 — c 2 sin 2 <p), where 

1-f Nsin 2 <p c 2 4-N 1 n-N 1 



v= 



; ■ zi i ■ ___: i ii : , 

(H-ftsin 2 <f>)(l — c 2 sin 2 <p) c 2 + w l-^-c 2 sin 2 <p • c 2 -f /& l + ftsin 2 ^ 



We can put %=1 and #=1 — <? 2 sin 2 <p, or vice versd^ as may be convenient. Making 
either substitution, and reducing to partial fractions the mean which we select for inte- 
gration, we have to multiply each fraction by the second value of % which doubles the 
number of fractions, and the index of their denominators ; but on again decomposing 
the fractions, and grouping them by their denominators, we reduce them to two more 
than those into which we had previously decomposed the mean. The two extra frac- 
tions will, of course, have 1+w sin 2 <p and 1 — c 2 sin 2 <p for their denominators. If we stop 
at the third arithmetic mean we shall thus have five partial fractions to integrate, and 
for the third harmonic mean, six. I do not actually exhibit the work or its results, 
because my doing so would not save labour to any one. Not only would the resulting for- 
mula be complicated with constants more easily managed in an arithmetical form, but it 
will seldom happen in practice that it will be worth while to reduce the elliptic function 
to the normal form given above. I have, however, done enough to show that my method 
is capable of approximately reducing any form, containing a function under the radical 
of the square root, to a small series of terms involving, at highest, logarithms or inverse 
tangents in their integrals. Moreover the approximation is so rapid, that, in the case 
of an elliptic integral of the third kind and of logarithmic form, nothing would be 
gained by having recourse to the interpolation of the only possible table, that of the 

double integral i — - j-~ = == ^ n \/l--c 2 sm 2 ().dQ 2 . The third pair of means will give 

six or more places of figures correct, and the fourth arithmetic mean is capable of giving 
twelve places. 

With respect to the actual integration of the partial fractions ultimately obtained, 
there is no difficulty. It will easily be seen that each partial fraction will be of the 

form z—. rr-. The integral of this with regard to <h is 

1 .- 

/T-— tan _i { v l-j~p.tan<p}. 

If 1+p is negative (=■— <f, suppose), this integral takes the form 

1 , fl 4-<7tan <p 

— 10£T v 

2c J be \l~~£tan<p 
MDCCCLX. 2 H 
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It is obvious that we may approximate in the same way to the values of ultra-elliptic 
integrals ; but the process will be more lengthy, on account of the greater complexity 
of these functions. 

In the case of higher radicals than the square root, so far as concerns a first approxi- 
mation only, it is clear that, if we insert several means between two integrable functions, 
any given geometric mean will be intermediate in value to the corresponding arithmetic 
and harmonic means ; but, inasmuch as the process affords no indication of what the 
second step is to be, it does not seem to have any useful application to such functions. 
But it brings all elliptic and ultra-elliptic functions within practical reach of the 
numerical computer. 
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for 



I have thought it advisable, upon reconsideration, to give the approximate formula 

v 1+Nsin 2 <p d(f> 



l+^sin 2 <p </].._ c 2 sin 2 <?5 



derived from the third mean of the arithmetic series :* 



IN fl.T 1 1 1 , I 



8 n 



8~8 n + c 2 ~ 4 w + £c an 4 2-f\/2 



"T A 



rc—N 



n-\- 



4 



n+. 



7^7T^ tan " 1{(1+ ^ 4tan?} 



•j 



1 N + c 2 1 



A A * I a 



+8-^T^7T=? tan " 1 ^ 1 - c ^- tan ^+IirH? 



1-^c 5 



-i 



tan -1 



-f 



1 



N + 



2-j-\/2 






2-h\/2 



— 1 



2 + s/2 «W 



4 



tan"" 1 



2 + \/2 



& 



tanp 



4 



+j 



jn -j- 



2 — \/2 



n + 



4 ^2-^ A 



■i 



tan 



-i 



2-v / 2 



tan?) 



i 



1— 2 c 



2 



tan<p 






1 also observe that com. log 

log 



2 + 



9-93123 06918 42 



2-^2 _ Q . 



=9-16567 93211 66. 



The application of this formula, in the shape given above, requires that 1+n be 
positive, and that <? sin p shall not exceed sin 45°, 
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If we equate n to N, we obtain for the approximate value of I (1 — & sin 2 <p)-%d<p, or F(c<p), 

Jo 



\ <p+\ (l-c 2 )-*tan- 1 {(l-c 2 )*tan?)}+j(l-^c 2 ) 'tan" 1 * 



i-y 



tan<p 



+i(l_£±v^ 

1 4 4 



"W 



► tan" 1 * 



+i 



4 



[(- 



2-f\/2 
4 

2~\/2 



tan<p 



c 2 ) tan<p 



l 

I now apply the numerical values p= -t, c= sin 45°, which, as has already appeared, 

lead to the most unfavourable case for approximation. For p=o*' a U the inverse 

1 
tangents become =oflr, and, by reducing the last two terms to a single radical, we 

easily obtain 

F(sin 45°, ~*r) =^{|+i sin 45° + J sin 60° + y^ 6 "^ 34 



34 



|tX 118034 09494 53=1-85407 52150. 



This exceeds the exact value given by Legendre, 1*85407 46773, by 0-00000 05377. 

The logarithm of the factor of \* is 0-07200 74743, of which the error is 0-00000 01290, 

in excess. 

It follows from this that the necessary error of the process, where the third mean of 
the arithmetic series is used, will never be more than a unit in the seventh place, and it 
is evident that this error will always be in excess. By an arbitrary correction, of which 
the amount may be easily guessed in any given case, the seventh place may always be 
made accurate. 



2 h 2 



